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ABSTRACT 


The  unique  properties  of  the  temperature  and  salinity  profiles  for  polar 
oeeans  are  eritical  for  high- latitude  mixed  layer  thermodynamics.  In  the  Polar  regions  the 
water  column  is  coldest  and  freshest  at  the  surface  where  ice  may  be  present.  This 
density  structure  often  leads  to  entrainment  and  affects  both  the  mixed  layer  depth  and 
the  ice  thickness. 

Thermobaricity,  the  combined  dependence  of  seawater  thermal  expansion 
on  temperature  and  pressure,  magnifies  the  buoyancy  flux  associated  with  mixed  layer 
convection.  When  thermobaricity  amplifies  entrainment  so  that  the  heat  into  the  mixed 
layer  is  greater  than  the  heat  leaving  the  water  column,  the  mixed  layer  warms  and  any 
existing  ice  begins  to  melt.  Similarly,  if  the  heat  entrained  is  less  than  the  heat  leaving 
the  column,  the  mixed  layer  cools  and  freezing  occurs  at  the  surface.  In  the  former 
situation  a  polynya,  or  region  of  no  ice  surrounded  by  ice  coverage,  may  form. 

A  one-dimensional  vertical  model  is  built,  and  trial  cases  are  run  to  show 
the  intricate  relationships  that  govern  the  heat  and  salt  fluxes  and  subsequent  ice 
thickness.  The  model  shows  the  importance  of  thermobaricity  to  the  air-sea-ice 
interactions.  It  also  offers  significant  insight  into  how  relatively  constant  atmospheric 
forcing  can  lead  to  polynya-like  conditions. 
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I.  INTRODUCTION 


Freezing  and  melting  of  iee  in  the  Aretie  and  Antaretie  environments  are 
determined  by  the  thermodynamie  fluxes  among  air,  sea-iee,  and  oeeanie  mixed  layer. 
The  energy  exehanged  among  these  three  domains  is  in  a  delieate  balanee  that  may 
theoretieally  lead  to  an  iee-free  Aretie  basin  or  another  iee  age  in  the  limits.  This  paper 
will  foeus  primarily  on  high-latitude  iee-mixed  layer  dynamics  in  both  the  northern  and 
southern  hemispheres.  Although  the  vertical  energy  fluxes  within  these  regions  are 
governed  by  similar  underlying  physics  one  dimensionally,  the  two  polar  regions  have 
differences.  One  phenomenon  possible  in  all  polar  regions  is  the  polynya. 

Polynyas  may  form  as  a  local  extreme  state  on  a  small  scale,  without  major  large- 
scale  or  basin  changes.  As  local  extremes,  polynyas  are  a  potential  laboratory  for  climate 
change.  Thus,  part  of  the  of  the  impetus  for  this  research  is  stated  clearly  by  Smith  et  al. 
(1990),  “At  predictable,  recurrent  locations  throughout  the  polar  regions,  there  are 
oceanic  areas  which  remain  either  partially  or  totally  ice  free  at  times  and  under 
climatological  conditions  where  we  would  expect  the  water  to  be  ice  covered.” 
Conversely,  ice  may  also  form  in  these  regions  at  times  when  climatological  conditions 
might  warrant  no-ice  conditions.  Simply  put,  there  is  more  to  ice  formation  and  melting 
than  climatological  conditions;  the  understanding  of  the  dynamics  and  thermodynamics 
of  the  water  column  below  and  coupled  with  the  ice  needs  to  be  explicitly  included  into 
prediction  models. 

This  work  has  many  implications  for  the  United  States  Navy  both  for  now  and  in 
the  future.  For  the  present  time  a  better  understanding  of  the  physics  of  the  mixed  layer 
will  improve  ice  thickness  and  ice  extent  forecast  models  like  PIPS  (Preller  et  al,  2002). 
The  output  from  such  a  model  could  be  vital  for  submarine  and  or  shipping  operations  in 
high-latitude  regions.  Furthermore,  information  gained  from  these  models  could  improve 
mission  effectiveness  and  efficiency  while  preventing  the  wasting  of  valuable  resources 
like  equipment,  time,  and  money.  For  example,  if  a  scenario  were  developed  where  a  P3- 
Orion  was  required  to  lay  sonobuoys  in  the  Barents  Sea  it  would  be  very  useful  in  the 
planning  stages  of  this  mission  to  have  information  regarding  the  ice  coverage.  If  the 
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region  is  ice  covered,  an  ice  penetrating  sonobuoy  would  be  very  useful  assuming  the 
thickness  of  the  ice  does  not  exceed  the  burying  capabilities  of  the  sonobuoy.  However, 
if  the  region  were  either  ice-free  or  the  ice  were  too  thick  then  the  employment  of  this 
type  of  sonobuoy  would  be  a  waste  of  time  and  equipment.  Similarly,  the  use  of  a 
conventional  sonobuoy  in  an  ice-covered  region  would  be  useless.  Either  way  real-time 
knowledge  about  the  ice  thickness  and  extent  would  be  invaluable  for  the  planning  of  this 
mission. 

The  long-term  importance  of  this  study,  for  the  Navy,  may  go  beyond  forecast 
capabilities  and  mission  effectiveness  to  include  macro-scale  circulation  and  deep  water 
formation.  The  same  physical  processes  that  are  explored  in  this  work  as  mechanisms  for 
ice  formation  and  melting  also  seemingly  play  a  significant  role  in  taking  oxygen  and 
nutrients  from  the  surface  and  reintroducing  them  back  into  the  water  column.  When 
these  nutrients  are  returned  deep  in  the  water  column  they  flow  with  the  denser  water  as  it 
spreads  around  the  globe  replenishing  the  world’s  oceans.  Understanding  completely 
how  this  process  works  and  the  resulting  ramifications  for  oceanic  and  atmospheric 
forecasts  could  provide  a  decisive  tactical  advantage. 

Chapter  II  explores  the  regional  descriptive  oceanography  of  the  high  latitudes. 
Climatology,  basin  circulation  patterns,  and  polynya  events  are  discussed  in  this  chapter. 
Following  this  background  discussion.  Chapter  III  delves  into  the  physics  of  the  coupled 
ice-mixed  layer  system,  examining  both  previous  scientific  studies  and  the  current  work. 
Requirements  for  the  accurate  modeling  of  this  physics  are  also  discussed.  Chapter  IV 
reviews  sensitivity  studies,  as  the  strengths  and  weaknesses  of  this  model  are  discovered. 
Finally,  Chapter  V  provides  a  summary  of  this  work  and  details  the  conclusions  drawn  as 
well  as  recommendation  for  work  that  needs  to  be  continued. 
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II.  POLAR  REGIONS 


This  section  is  designed  to  provide  the  background  knowledge  that  is  important 
for  the  air-sea-ice  model.  The  focus  of  the  circulation  section  is  on  the  Arctic  since  the 
flow  patterns  in  the  Northern  Hemisphere  are  much  more  complex  than  the  circulation 
patterns  in  the  Southern  Hemisphere.  It  is,  however,  important  to  note  that  the  Weddell 
Polynya  in  Antarctica  appeared  to  migrate  west  from  year  to  year  during  its  presence  in 
the  1970s.  Despite  this  fact,  for  the  purposes  of  the  one-dimensional  vertical  model,  it  is 
assumed  that  if  the  model  can  handle  the  horizontal  advective  effects  of  the  Arctic,  it  is 
robust  enough  to  handle  similar  effects  in  the  Antarctic.  In  regards  to  climate,  similar 
assumptions  cannot  be  made  since  the  atmospheric  forcing  can  be  significantly  more 
extreme  in  the  Antarctic  than  in  the  Arctic.  The  section  concerning  leads  and  polynyas 
focuses  primarily  on  the  Weddell  Polynya  in  the  Antarctic.  The  reason  for  this  is  not 
because  physics  does  not  support  polynyas  in  the  Arctic,  but  simply  because  the 
relatively  large  amount  of  information  about  the  Weddell  Polynya  makes  it  a  good 
example  of  and  laboratory  for  the  physics  that  create  ice-free  regimes  in  an  ice-covered 
region. 

A,  CIRCULATION  OF  THE  ARCTIC  OCEAN 

To  better  understand  one-dimensional  vertical  effects  in  the  Arctic  Ocean,  the 
horizontal  movement  of  water  into  and  through  the  region,  as  well  as  the  bathymetry  of 
the  region  need  to  be  discussed  (Fig.  1).  For  the  purposes  of  this  paper  the  rather 
complex  circulation  patterns  and  bottom  topography  of  the  Arctic  Ocean  examined  by 
Rudels  et  al.  (1994)  is  simplified  to  include  two  primary  basins,  the  Eurasian  Basin  and 
the  Canadian  Basin,  the  Lomonosov  Ridge,  and  the  flow  of  water  between  and  around 
these  features.  This  simplification  is  justified  since  Rudels  et  al.  (1994)  demonstrate  a 
strong  front  between  the  two  basins  at  the  Lomonosov  Ridge.  Similarly,  Aagaard  and 
Carmack  (1994)  show  evidence  to  support  the  presence  of  a  strong  boundary  current 
around  the  Lomonosov  Ridge. 
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Figure  1.  The  Aretie  Oeean  with  bathymetry.  From  Barry  et  al,  1993. 
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Although  water  can  enter  and  leave  the  Arctic  Ocean  from  several  points,  the 
Bering  Strait  will  be  considered  the  primary  point  of  entry  for  Pacific  water  per  Aagaard 
and  Carmack  (1994).  Likewise,  the  Fram  Strait  was  considered  the  main  source  for 
Atlantic  and  Intermediate  waters  according  to  Rudels  et  al.  (1994).  A  secondary  source 
of  water  into  the  Arctic  is  river  runoff  that  can  act  to  cap  a  water  column  and  prevent  heat 
exchange  with  the  environment  (Mellor  and  Hakkinen,  1994). 

Despite  their  proximity  to  one  another,  the  two  basins  are  different.  Rudels  et  al. 
(1994)  suggest  that  the  circulation  through  the  Canadian  Basin  is  slower  than  that  of  the 
Eurasian  Basin,  causing  the  water  to  reside  in  the  Canadian  Basin  by  as  much  as  a  decade 
longer  than  in  the  Eurasian  Basin.  Rudels  et  al.  (1994)  show  the  Canadian  Basin  to  be 
well  mixed  according  to  temperature  and  salinity  by  water  type.  Rudels  et  al.  (1994)  also 
depict  the  Eurasian  Basin  as  being  well  mixed  according  to  temperature  and  salinity  by 
water  type.  However,  the  Atlantic  and  intermediate  waters  of  the  Eurasian  Basin  increase 
in  both  temperature  and  salinity  as  depth  increases,  whereas  the  same  type  of  water  in  the 
Canadian  Basin  decreases  in  both  temperature  and  salinity  as  depth  increases  (Rudels  et 
al.,  1994).  This  difference  is  important  to  the  vertical  heat  flux  and  convection  of  each 
basin  respectively.  Similarly,  the  difference  is  important  when  water  from  both  basins 
come  together  at  the  Eomonosov  Ridge  because  double  diffusion  may  occur  (Rudels  et 
al.,  1994). 

The  circulation  time  is  shorter  in  the  Eurasian  Basin  than  the  Canadian  Basin 
because  advection  is  strong  enough  in  the  Eurasian  Basin  to  pass  the  incoming  water 
back  out  of  the  Eram  Strait  (Rudels  et  al.,  1994).  That  does  not  mean  that  vertical 
processes  do  not  occur  in  this  basin.  In  fact,  the  temperature  and  salinity  profiles  detailed 
by  Rudels  et  al.  (1994)  could  allow  entrainment  of  warmer  underlying  water  into  the 
colder  surface  layer.  If  this  were  to  occur,  the  rising  heat  could  melt  ice.  At  the  same 
time,  the  loss  of  heat  at  depth  would  cool  the  lower  layers. 

The  temperature  and  salinity  structure  of  the  Canadian  Basin  is  of  particular 
interest  since  Rudels  et  al.  (1994)  suggest  that  the  structure  will  allow  heat  to  be  mixed 
down  into  the  water  column  (Eig.  2).  Mixing  of  heat  down  into  the  water  column  would 
deepen  the  mixed  layer. 
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Figure  2.  Representative  temperature  and  salinity  profiles  from  various  loeations  in  the 
Northern  Hemisphere  high  latitudes.  From  Coaehman  and  Aagaard,  1974. 
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If  this  mixing  of  heat  down  is  greater  than  the  entrainment  up  from  the  abyss,*  or 
if  the  basin  is  not  deep  enough  for  an  abyss  to  exist,  then  little  heat  would  be  added  to  the 
mixed  layer  to  melt  ice.  If,  however,  ice  forms  because  entrained  heat  is  less  than  heat 
lost  to  the  atmosphere,  the  mixed  layer  salinity  increases.  Increased  salinity  would  then 
make  the  water  sink  convectively.  This  process  creates  the  warm  and  salty  water  at  depth 
discussed  by  Aagaard  and  Carmack  (1994). 

Introducing  advective  processes  like  river  runoff  or  narrow  currents  like  the  West 
Spitsbergen  Current,  as  mentioned  by  Mellor  and  Hakkinen  (1994),  probably  do  little  to 
offset  the  vertical  processes  in  the  Canadian  Basin.  River  runoff  would  create  fresh  water 
at  the  surface  that  would  quickly  freeze  and  add  to  the  ice.  Currents  at  either  the  surface 
or  some  depth  could  create  eddies  and  cause  stirring  in  the  horizontal.  However,  due  to 
the  stronger  mixing  of  heat  down  in  the  column,  these  eddies  would  have  little  effect  on 
the  direction  of  vertical  mixing. 

Advective  processes  like  friction,  wind  stress,  and  currents  undoubtedly  keep  the 
ice  in  the  Canadian  Basin  from  growing  unbounded.  Similarly,  Lemke  et  al.  (1990)  point 
out  that  a  deep  mixed  layer  can  resist  ice  growth.  This  situation  would  be  very  similar  to 
the  circumstances  surrounding  the  presence  of  a  polynya.  However,  strong  mixing  of 
heat  down  in  the  Canadian  Basin  could  theoretically  make  it  a  one-dimensional  ice 
generator,  especially  when  the  temperature  of  the  water  at  the  surface  is  at  or  near  the 
freezing  point.  As  it  is,  the  Canadian  Basin  becomes  a  density  capacitor,  building  up 
density  over  longer  periods  of  time  and  then  slowly  but  constantly  discharging  salty  and 
relatively  warm  water  at  the  Lomonosov  Ridge.  This  water  then  mixes  with  cooler  and 
fresher  water  released  from  the  Eurasian  Basin.  When  the  waters  from  the  two  basins 
come  in  contact  at  the  Lomonosov  Ridge  double  diffusion  may  occur  as  noted  above. 
Strong  boundary  currents  along  the  Lomonosov  Ridge  and  the  east  coast  of  Greenland 
then  push  and  carry  the  water  before  it  exits  from  the  Lram  Strait  at  depth  (Aagaard  and 
Carmack,  1994).  Stable  water,  defined  as  cold  salty  water  beneath  warm  fresher  water,  is 
the  ultimate  result  of  the  double  diffusion.  In  this  way  warm,  salty,  cold,  and  fresh  waters 


*  The  term  abyss  is  borrowed  from  Hakkinen  (1987)  and  is  used  here  to  mean  the  depth  of  the  water 
column  beneath  the  mixed  layer. 
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are  all  brought  into  the  Aretie  Oeean,  split  and  proeessed  by  both  the  topography  and 
eireulation,  and  then  re-circulated  to  replenish  the  world’s  oceans  (Fig.  3). 


B,  ATMOSPHERIC  CONDITIONS 

Atmospheric  forcing  plays  a  role  in  both  the  northern  and  southern  hemispheres. 
Ultimately,  ice  will  form  or  melt  as  a  result  of  the  overall  energy  balance.  This  balance 
includes  radiation  fluxes,  sensible  heat,  and  latent  heat  (Fig.  4)  in  both  polar  regions. 
Whether  open  areas  like  leads  exist,  when  heat  is  brought  to  the  surface,  is  critical  to  the 
energy  balance  since  Vihma  (1995)  states  that  ice  cover  prevents  heat  exchange.  Okada 
and  Yamanouchi  (1995)  similarly  show  that  when  the  sea  ice  concentration  increases  in 
the  Antarctic  winter  the  surface  heat  flux  decreases.  If  a  large  section  of  ice  in  the  South 
Pole  were  instantaneously  removed,  evidence  presented  by  Wu  et  al,  (1996)  suggests 
that  the  energy  balance  would  restore  the  ice  within  one  year.  This  not  only  demonstrates 
the  strength  of  the  energy  balance  but  also  suggests  a  seasonality  to  ice  formation. 

Seasonal  change  is  also  an  important  factor  in  the  overall  energy  balance.  Sea  ice 
concentration  is  greatest  in  the  southern  hemisphere  in  the  July  to  August  time  frame 
during  the  Antarctic  winter  (Fig.  5).  Similarly,  the  sea  ice  extent  is  greatest  in  the 
northern  hemisphere  from  February  to  March  during  the  Arctic  winter  (Fig.  6).  Seasonal 
differences  can  be  witnessed,  and  melting  and  freezing  can  be  observed  through  surface 
albedo.  A  low  albedo,  measured  in  percentage,  corresponds  to  melting  ice  and  a  high 
albedo  classifies  solid  ice  (Barry  et  al,  1993). 

There  are  also  significant  seasonal  differences  in  ice  coverage  within  each  of  the 
polar  regions.  Preller  et  al  (2002)  describe  ice  coverage  in  Arctic  regions  to  be  thinnest 
during  the  summer;  from  June  to  September,  and  thickest  in  the  winter;  from  March  to 
April.  Similarly,  Preller  et  al  (2002)  suggest  that  different  seas,  within  the  Arctic,  for  a 
given  season  will  have  different  ice  thickness  and  coverage,  but  that  no  sea  will 
perpetually  have  the  greatest  coverage.  Preller  et  al  (2002)  compare  this  phenomenon  to 
a  seesaw  where  for  a  given  summer  sea  “A”  has  thicker  ice  than  sea  “B”,  but  the 
following  summer  sea  “B”  will  have  greater  thickness  than  sea  “A”. 
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Figure  3.  Circulation  of  the  Atlantic  Layer  and  intermediate  depth  waters  of  the  Arctic 

Ocean.  From  Rudels  et  a/.,  1994. 
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SURFACE  ENERGY  BUDGET  (60S-70S) 

(W/rrf) 


Figure  4.  Typical  annual  surface  energy  budget  for  the  Southern  Hemisphere.  From  Okada 

and  Yamanouchi,  1995. 
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Figure  5.  Typical  annual  Southern  Hemisphere  sea  ice  concentration.  From  Okada  and 

Yamanouchi,  1995. 
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Figure  6. 


Typical  annual  Northern  Hemisphere  sea  ice  extent.  From  Walsh  and  Johnson, 

1979. 
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Understanding  both  seasonal  and  interannual  seesaw  effeets  are  important  for  the  United 
States  Navy  when  attempting  to  foreeast  iee  thiekness  and  extent.  To  do  this,  the  Navy 
uses  the  Polar  lee  Predietion  System  (PIPS),  (Preller  et  al,  2002).  This  model  makes  use 
of  real  time  and  elimatologieal  atmospherie  data,  the  Cox  Oeean  Model,  and  the  Hibler 
lee  model  to  produee  a  foreeast  for  iee  approximately  five  days  in  the  future  (Preller  et 
al,  2002)  (Fig. 7).  One  problem  with  PIPS  doeumented  by  Preller  et  al  (2002)  is  that  it 
often  over-prediets  iee  in  the  Barents  Sea  and  under-prediets  iee  in  the  Labrador  Sea. 

This  over-predietion  of  iee  may  be  attributed  to  the  poor  resolution  of  smaller  adveetive 
proeesses  like  eurrents  or  river  run-off  (Mellor  and  Hakkinen,  1994).  Another  possibility 
for  the  diserepaneies  in  predietion  is  that  the  model  does  not  resolve  the  heat  fluxes  in 
these  regions  as  aeeurately  as  needed.  Speeifieally,  the  Labrador  Sea  has  heat  fluxes  and 
deep  mixed  layers  more  eonsistent  with  the  Antaretie  region  than  other  Aretie  seas,  and 
so  the  vertieal  proeesses  of  this  sea  generate  greater  amounts  of  iee  than  predieted  (Guest, 
personal  oommunieation). 


C.  LEADS  AND  POLYNYAS 

Leads  and  polynyas  are  both  areas  of  open  water  where  iee  would  have  been 
predieted  given  the  elimatology  of  the  surfaee  foreing  by  the  atmosphere.  There  are  two 
types  of  polynyas  enumerated  by  Smith  et  al  (1990);  the  first  is  a  “latent  heat”  polynya 
that  is  maintained  by  horizontal  adveetive  heating  proeesses  attributable  to  wind  and 
eurrents.  The  seeond  is  a  “sensible  heat”  polynya  that  is  ereated  from  vertieal  eonveetive 
heating  from  the  oeean.  In  either  ease  the  net  heat  flux  into  that  area  is  suffieient  to 
inhibit  iee  formation.  As  this  researeh  is  foeused  primarily  on  vertieal  effeets,  the 
“sensible  heat”  polynya  is  of  greatest  interest.  Smith  et  al  (1990)  deseribe  leads  as  open 
ehannels  within  the  iee  that  are  ereated  by  a  divergenee  in  the  iee  flow.  Leads  are  also  of 
interest  espeeially  when  eonsidering  that  oeean  dynamies  and  heating  from  below  may 
develop  weaknesses  in  the  iee.  When  the  weaker  iee  is  put  under  stress  from  wind  or 
surrounding  iee  it  diverges  into  open  water  ehannels. 
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Polar  Ice  Prediction  System  2.0 


Figure  7.  U.S.  Navy  PIPS  Model.  From  Preller  et  al, 
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2002. 
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While  lead  location  does  have  a  stochastic  nature,  polynyas  often  recur  in  the 
same  location  with  some  periodicity.  However,  it  is  important  to  note,  as  Smith  et  al. 
(1990)  does,  that  the  Weddell  Sea  sensible  heat  polynya  disappeared  in  the  1980’s  after 
having  appeared  regularly  in  the  late  1970’s.  This  indicates  the  strong  dependence  of 
open  water  formation  on  the  dynamics  occurring  beneath  the  surface  and  suggests  that 
climatological  conditions,  while  important,  are  not  the  sole  causes.  Similarly,  the 
seasonal  nature  of  the  heat  flux,  with  greater  heat  flux  out  of  the  water  in  winter,  is  no 
assurance  of  ice  formation;  nor  does  summer  warming  of  the  water  column  guarantee 
complete  ice  melting. 

In  regards  to  the  Weddell  Sea  in  Antarctica,  Martinson  et  al.  (1981)  assert  that  the 
Weddell  Polynya  cannot  be  created  by  wind  forcing  or  atmospheric  effects  alone  since 
the  low-pressure  system  required  to  create  it  would  not  survive  long  enough  in  this 
environment.  Instead,  they  postulate  that  some  type  of  preconditioning  occurs  to  raise 
the  pycnocline  from  below  the  mixed  layer.  Martinson  et  al.  (1981)  offer  several 
possibilities  for  the  source  of  the  preconditioning  including  both  topography  and  mean 
flow  spin-up.  The  preconditioning,  whatever  its  source,  coupled  with  salt  rejection  at  the 
surface  during  freezing  causes  a  layer  of  cold  fresh  water  above  a  raised  portion  of  warm 
salty  water  from  below.  When  the  density  of  the  cold  fresh  water  surpasses  that  of  the 
warm  salty  water  the  column  can  overturn  leaving  warm  salty  water  at  the  surface.  In 
this  way  freezing  of  ice  initiates  a  negative  feedback  where  the  end  result  is  ice  melting 
from  the  heat  of  the  newly  raised  water  just  below  the  ice  at  the  surface  of  the  water 
column.  Because  of  density  driven  overturning  Martinson  et  al.  (1981)  contend  that  the 
temperature  and  salinity  profiles  of  both  the  pycnocline  beneath  the  mixed  layer  and  the 
mixed  layer  itself  are  vital  for  determining  the  formation  of  the  Weddell  Polynya.  It  is  of 
interest  to  note,  as  well,  that  although  no  conclusions  are  made  in  regards  to  the  effects  of 
topography,  Martinson  et  al.  (1981)  do  speculate  that  the  polynya  is  subject  to  the 
circulation  of  the  Weddell  gyre  and  appears  to  cycle  from  the  east  to  the  west  with  the 
mean  flow.  Furthermore,  they  suggest  evidence  that  the  polynya  moves  to  an  extreme 
westward  location  before  it  is  recreated  at  its  extreme  eastern-most  point  (Martinson  et 
al,  1981). 
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III.  THEORY 


A.  PREVIOUS  STUDIES 

Recognizing  the  importance  of  atmospheric  conditions  and  seasonal  variations  is 
necessary  for  predicting  ice  freezing  and  melting.  Similarly,  as  in  the  case  of  the 
Weddell  Polynya  and  mentioned  by  Martinson  et  al.  (1981),  the  circulation  patterns,  both 
surface  and  subsurface,  and/or  bathymetry  may  aid  or  precondition  melting  or  freezing 
effects.  However,  since  low  or  high  pressure  systems  and  their  accompanying  winds 
come  and  go  at  varying  time  scales,  and  bathymetry  remains  relatively  constant  within  a 
region,  the  physics  of  the  water  column  is  the  most  significant  component  of  sustained 
surface  freezing  and  melting.  More  specifically,  the  coupled  thermodynamics  and 
dynamics  that  together  mix  water  vertically  upwards  and  downwards  contribute  the  most 
to  ice  freezing  or  melting.  Lemke  et  al.  (1990)  explain  the  importance  of  this  coupling  by 
suggesting  that  a  situation  may  develop  where  ice  melts  due  to  thermodynamic  effects 
creating  an  area  of  no  ice,  but  then  freezes  again,  almost  immediately,  due  to  water 
movement  at  the  surface.  This  argument  highlights  the  need  for  the  source  of  the  heat,  in 
addition  to  the  heat  it  provides,  to  be  brought  to  the  surface  for  a  sustained  no-ice  area  to 
exist  in  the  polar  climate. 

Hakkinen  (1987)  explains  the  thermodynamics  and  dynamics  necessary  to  bring 
heat  to  the  surface  and  sustain  it  using  a  two-layer  model.  The  two  layers  suggested  by 
Hakkinen  (1987)  are  a  cold  fresh  surface  layer  and  a  warm  salty  bottom  layer  that  is  at 
rest.  As  described  by  Kraus  and  Turner  (1967)  the  heat  and  salinity  fluxes  between  these 
two  layers  are  dependent  upon  the  difference  between  the  temperature  and  salinity 
profiles  of  each  layer  respectively.  Hakkinen  (1987)  investigates  the  vertical  transfer  of 
heat  on  the  mesoscale.  Specifically,  she  focuses  on  entrainment  of  warm  salty  water 
from  the  lower  layer  up  into  the  surface  layer  at  the  ice  edge  in  the  marginal  ice  zone 
(Hakkinen,  1987).  Her  one-dimensional  model  suggests  that  entrainment  has  an  inverse 
dependence  on  density;  as  the  depth  of  the  mixed  layer  increases  entrainment  decreases, 
thereby  ensuring  that  the  mixed  layer  does  not  deepen  infinitely  (Hakkinen,  1987). 
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Dynamically,  Hakkinen  (1987)  relies  on  strong  winds  blown  aeross  the  surfaee  to 
ereate  drag  on  the  water  eolumn.  This  drag,  whieh  is  three  times  as  great  when  there  is 
iee  at  the  surfaee  as  when  there  is  no  ice  present,  creates  strong  Ekman  transport 
(Hakkinen,  1987).  Not  surprisingly,  the  greater  the  drag  the  greater  the  Ekman  transport, 
and  so  Ekman  transport  is  greatest  under  the  iee  (Hakkinen,  1987).  In  order  to  quantify 
the  drag  Hakkinen  (1987)  defines  a  fraetional  ealeulation  of  iee  eompaetness,  A. 
Hakkinen’ s  (1987)  preeonditioning,  as  termed  by  Martinson  et  al.  (1981)  for  the  priming 
of  upeoming  interaetion  between  the  two  layers,  eomes  from  the  stronger  Ekman 
transport  down  beneath  the  iee  that  leads  to  inereased  upwelling  at  the  iee  edge.  Through 
this  meehanism  upwelling  assists  entrainment  in  keeping  warmer  and  more  saline  water 
at  the  surface  and  keeps  the  area  free  of  iee  (Hakkinen,  1987). 

Eemke  et  al.  (1990)  similarly  use  a  one-dimensional  two-layer  model  eonsistent 
with  Kraus  and  Turner  (1967)  to  explain  sustained  heat  at  the  surfaee.  Eike  Hakkinen, 
Eemke  et  al.  (1990)  is  eoncerned  with  the  momentum  balance  between  the  atmosphere, 
iee,  and  water  column.  Eemke  et  al.  (1990)  eombines  the  equations  of  motion,  iee 
eompaetness  and  eoverage,  the  nonlinear  viseous-plate  eonstitutive  law  from  Hibler 
(1979),  and  atmospherie  foreing  to  explain  vertical  movement  of  heat.  Eemke  et  al. 
(1990)  also  find  that  the  greatest  vertieal  heat  flux  in  the  Weddell  Sea  is  at  the  iee  edge. 

Unlike  Hakkinen  (1987),  however,  Eemke  et  al.  (1990)  focus  less  on  atmospherie 
forees  like  wind  and  more  on  the  entrainment  veloeity  as  the  mechanism  for  heating  at 
the  surface.  In  other  words,  the  overall  oeeanic  heat  flux  is  dependent  on  the  entrainment 
veloeity  (Eemke  et  al.,  1990).  This  means  that  whether  ice  melts  or  freezes  it  is  directly 
dependent  on  the  entrainment  veloeity.  Aeeording  to  Eemke  et  al.  (1990)  positive 
entrainment  brings  heat  to  the  surfaee  and  melts  the  iee,  entrainment  less  than  zero 
extends  the  mixed  layer  to  the  depth  of  the  water  eolumn,  and  entrainment  equal  to  zero, 
setting  the  mixed  layer  depth  equal  to  the  Mon  in -Obukhov  length,  oecurs  when  the  two 
layers  are  thermodynamieally  independent  of  one  another.  Eemke  et  al.  (1990)  formulate 
the  entrainment  values  to  be  dependent  upon  the  equation  of  state,  and  thereby  effectively 
emphasize  the  role  of  the  temperature  and  salinity  expansion  eoefficients.  Eurthermore, 
the  inelusion  of  the  equation  of  state  underscores  the  importance  of  ehanges  in  both 
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temperature  and  salinity  between  the  abyss  and  the  mixed  layer  for  vertieal  heat  exehange 
(Lemke  et  al,  1990). 


B,  THERMOBARICITY 

Thermobarieity  is  an  important  proeess  to  inelude  when  attempting  to  model  or 
understand  the  eomplete  physies  of  the  mixed  layer.  It  is  defined  as  the  eombined 
dependenee  of  seawater  thermal  expansion  on  temperature  and  pressure  (MeDougall, 
1984).  The  work  done  by  Lemke  et  al.  (1990)  lays  the  groundwork  for  the  signifieanee 
of  thermobarieity.  It  is  important  beeause  it  magnifies  the  buoyaney  flux  assoeiated  with 
mixed  layer  eonveetion  (Garwood,  1991).  If  thermobarieity  is  ineluded  in  the  equation  of 
state  the  thermal  expansion  eoeffieient,  a  ,  is  no  longer  a  eonstant,  but  is  instead  a 
funetion  of  both  temperature  and  depth.  This  means  that  entrainment  is  not  only  a 
funetion  of  the  differenees  in  temperature  and  salinity  between  the  abyss  and  mixed  layer, 
but  is  also  a  funetion  of  the  depth  and  temperature  of  the  abyss-mixed  layer  boundary. 

The  physieal  result  of  the  inelusion  of  thermobarieity  is  eonsistent  with  Lemke  et  al. 
(1990);  when  entrained  heat  from  below  is  greater  than  the  heat  leaving  the  water 
eolumn,  the  mixed  layer  warms  and  any  existing  iee  begins  to  melt.  Similarly,  if 
entrained  heat  is  less  than  the  heat  leaving  the  eolumn  the  mixed  layer  eools  and  freezing 
oeeurs  at  the  surfaee.  However,  the  magnitude  of  the  entrainment  eould  be  very  different 
with  the  inelusion  of  thermobarieity,  under  certain  conditions,  resulting  in  either  freezing 
or  melting  when  Lemke  et  al.  (1990)  predict  the  opposite. 

Garwood  et  al.  (1994)  shows  the  necessity  of  understanding  thermobarieity  and 
the  importance  of  being  able  to  model  it.  He  enumerates  four  different  scenarios  where 
thermobaric  effects  have  implications  in  polar  modeling. 

1)  Mixed-layer  turbulent  kinetic  energy  (TKE)  may  be  increased  several-fold  by 
the  thermobaric  amplification  of  the  buoyancy  flux  (Garwood,  1991).  This  phenomenon 
is  not  significant  at  temperate  and  low  latitudes,  but  it  becomes  large  for  polar  sea  mixed 
layers.  It  cannot  be  neglected  for  realistic  air-sea  ice  interactions.  Large  eddy  simulation 
(LES)  verifies  its  importance  in  increasing  entrainment.  However,  commonly  used 
parameterizations  of  ocean  mixing,  including  second-order  closure  of  Mellor  and 
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Yamada  (1982)  and  k-profile  parameterization  (KPP)  (Large  et  al.,  1994),  do  not  yet 
include  this  physics. 

2)  Buoyancy  is  not  conserved  in  coastal  subduction.  Although  buoyancy  is  well 
conserved  in  vertical  advection  of  water  masses  at  lower  latitudes,  it  is  not  a  conservative 
property  in  the  polar  seas.  This  is  especially  important  for  vertical  circulation  that  is 
enhanced  on  shelf  slopes.  Here  it  becomes  a  severe  limitation  for  layered  models  that 
presume  a  constant  potential  density  for  each  layer  (Garwood  and  Isakari,  1993;  Jiang 
and  Garwood,  1995,  1996,  1998). 

3)  "Thermobaric  stability"  is  an  increase  in  the  hydrostatic  stability  of  seawater 
having  positive  spiciness  and  nearly  density-compensating  T-S  structure.  It  causes  spicy 
warmer  and  saltier  parcels  to  resist  vertical  mixing  because  of  a  newly  recognized 
thermobaric  component  to  the  buoyancy  frequency.  This  added  component  depends 
upon  the  thermobaric  depth  .  This  phenomenon  enables  internal  gravity  waves 

without  vertical  buoyancy  gradient  (Garwood,  2002),  and  it  helps  explain  lateral  intrusion 
of  filaments  of  spicy  water  across  T-S  frontal  zones.  Furthermore,  the  theory  predicts 
reduced  stability  for  cold  and  fresh  filaments  with  negative  spiciness.  This  process  is 
potentially  important  in  enhancing  shoreward  transport  but  reducing  the  flux  of 
suspended  and  dissolved  material  from  the  shelf  into  the  Arctic  interior. 

4)  A  critical  mixed-layer  depth  for  maintenance  of  open  water  and  sensible 

heat  polynyas  (Smith  et  al,  1990)  under  freezing  atmospheric  temperatures  is  predicted 
because  of  thermobaricity  (Garwood,  2001).  Whenever  the  mixed  layer  depth  exceeds  a 
critical  depth,  the  entrainment  heat  flux  will  exceed  the  surface  heat  loss  to  the 

atmosphere  ,  resulting  in  ice  melting  regardless  of  the  rate  of  heat  loss  to  the 

atmosphere  from  the  ice-mixed  layer  system.  This  is  a  positive  feedback  mechanism  for 
warming  the  surface  layer,  limiting  ice  formation,  and  cooling  the  ocean  interior.  It 
undoubtedly  plays  a  role  in  establishing  the  extent  and  geographic  evolution  of  the 
marginal  ice  zone,  and  in  the  maintenance  of  polynyas  and  their  entrainment  of  nutrient- 
rich  water  into  the  surface  layer. 
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C.  FUNDAMENTAL  EQUATIONS 

Since  Kraus  and  Turner  (1967)  first  modeled  the  mixed  layer,  oeeanographers 
have  worked  to  understand  the  role  of  turbulenee  in  mixed  layer  dynamies.  Their  foeus 
was  limited  to  the  temperate  oeeans  and  the  role  of  air-sea  interaetions  on  turbulenee. 
Independently,  sinee  Hibler  (1979),  polar  oeeanographers  have  proposed  and  advanced 
models  for  iee  mechanics,  with  cursory  attention  to  mixed  layer  turbulenee.  Recognizing 
the  need  to  eonsider  mixed  layer  dynamies  under  the  ice,  Lemke  et  al.  (1990)  and 
Hakkinen  (1987)  included  linear  formulations  for  entrainment  into  their  coupled  iee- 
oeean  numerieal  models.  However,  Garwood  (1991)  showed  that  the  equation  of  state’s 
nonlinear  properties  needed  to  be  ineluded  when  eonsidering  turbulent  mixing  in  the  eold 
water  of  the  polar  seas. 

The  physieal  proeesses  of  the  eoupled  system  depend  heavily  on  the  vertieal 
entrainment  of  warmer,  saltier,  and  denser  water  from  below.  Normally  in  the  lower 
latitudes  the  denser  water  from  below  is  both  eolder  and  saltier  than  the  water  above  it, 
but  in  the  high  latitudes  where  salinity  is  the  dominate  determinate  of  density,  it  is  the 
warm  and  saline  water  that  makes  up  the  high  density  water  from  below.  Entrainment,  or 
the  rate  of  deepening  of  the  mixed  layer,  is  important  for  two  primary  reasons.  First, 
entrainment  is  needed  to  ealculate  both  the  temperature  and  salinity  fluxes  that  make  up 
the  buoyaney  flux.  Seeond,  when  ice  is  present  at  the  surfaee,  entrained  water  may  heat 
the  ice  layer  directly  without  heating  the  mixed  layer  significantly,  and  so  ice  thickness  is 
directly  dependent  on  the  entrainment. 


Prediction  of  polar  sea  mixed  layer  dynamies  and  thermodynamies  is  dependent 
upon  the  same  physies  that  governs  all  water  pareel  motion:  the  Reynolds  equations  for 
the  mean  momentum,  which  are  derived  from  the  Navier-Stokes  equations, 

U^+UU^+Wy+WU^=-  —  +  fv  -  fyW-{u'u''^  -{u'v'^  -{u'w'^  +v{uxx+Uyy+Uzz^{\) 


u,+uu+vu^,+  wu^  = - -  + 

y  ^ 


w,  -I-  uw^  -I-  wv„  -I-  ww,  = - ^  -I- 


J^-fyW-(u'u')  -(u'v')  -(u'w')  +v( 


Kxx  H  yy  Hzz 


)(2) 


T 


P 


fyU-g-{u'w'^  -(vV)  -(wV)  +v{wxx+Wyy  +Wzz'^{2)) 


the  continuity  equation. 
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(4) 


u,+Vy+w^=0 

the  heat  and  salt  equations  where  the  Reynolds  deeomposition  is  also  applied  with 
T  =  T  +  T',S  =  S  +  S' ,  averaged,  and  put  into  flux  form  as, 

!^  =  T,  +  uT,  +  W^+wT^  =-(T^\  -  (^),  -  (^'),  +  *'r  (t  +  +  ?; )  (5) 

=  S,  +ii5,  +  W,  +  iv J  =  -(sv)^ - (SV')^  - (SV')^  +  *•, (s„  +  +  4 )  (6) 

and  the  equation  of  state, 

p  =  pM-a(T-Tj  +  Ji(S-Sj]  (7) 


Predietion  of  the  turbulent  fluxes,  and  entrainment  in  partieular,  is  based  upon  the 
Turbulent  Kinetie  Energy  (TKE)  equation,  whieh  is  also  derived  from  the  Navier-Stokes 
equations. 
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This  equation  ealeulates  the  energy  available  to  do  work  ereated  from  shear  produetion, 

,  and  the  vertieal  fluxes  of  temperature  and  salinity  that  make  up  the 


,du 

dz 


U'w' - h  V  V - 


dz 


vertieal  buoyaney  flux,  b'w'  =  agT'w'-  /3gS'w' .  The  four  terms  to  the  right  of  the 
equals  sign  in  the  TKE  equation  are,  from  left  to  right,  as  follows;  the  shear  produetion  of 
horizontal  turbulent  kinetie  energy,  the  production  or  damping  of  vertical  turbulent 
kinetic  energy  by  buoyancy  forces,  redistribution,  and  dissipation.  These  terms  apply  to 
the  mixed  layer  whether  or  not  ice  is  present  at  the  surface.  The  shear  production  term  is 
always  non-negative,  and  horizontal  turbulent  kinetic  energy  is  a  direct  result  of  the  wind 
stress  at  the  surface  and  produces  “forced  convection.”  Vertical  turbulent  kinetic  energy 
is  produced  whenever  the  buoyancy  flux  is  positive  due  to  negative  salinity  flux  and/or 
positive  temperature  flux.  Positive  buoyancy  flux  results  in  “free  convection.” 

Buoyancy  flux  will  always  be  negative  at  the  bottom  of  the  mixed  layer,  where  the 
turbulent  kinetic  energy  is  used  to  entrain  the  underlying  heavier  water  into  the  mixed 
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layer.  Redistribution  moves  turbulent  kinetie  energy  vertically  within  the  mixed  layer. 
Dissipation  converts  turbulent  kinetic  energy  to  heat,  and  is  always  a  loss  term. 


To  solve  for  the  entrainment  it  is  necessary  to  solve  for  the  vertical  integral  of  the 
TKE  equation.  If  equation  (8)  is  vertically  integrated  over  the  mixed  layer,  the 
entrainment  velocity  may  be  calculated  in  terms  of  the  stratification  (AT  and  A^),  mixed 
layer  depth,  h,  and  the  forcing  velocity  scales  (w,  and  w, ),  which  are  defined  below. 
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The  resulting  equation  (9)  shows  that  all  of  the  available  energy  to  do  work,  or  stirring  in 
this  case,  is  provided  from  either  forced  convection  from  wind  or  ice  forcing,  or  from  the 
free  convection  from  within  the  water  column.  Forced  convection. 


CjM, 


=  1 


—r-,du  —r-,dv 
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(10) 


originates  in  the  atmosphere  and  quantifies  the  turbulence  created  from  wind  blowing 
across  the  surface  or  across  ice  at  the  surface.  Free  convection,  or  positive  c^wi' ,  is  the 
vertical  integral  of  the  buoyancy  flux  over  the  mixed  layer, 

=  f  b'w'dz  (11) 

J  -h 

As  can  be  seen  in  the  equation  for  entrainment  (9)  the  equation  of  state  (7)  has  a 
significant  influence  on  the  buoyancy  flux  and  magnitude  of  free  convection  (11).  It 
should  be  mentioned  that  the  term  “free  convection”  only  applies  when  wj’  is  positive.  If 

is  negative  the  integrated  buoyancy  flux  dampens  the  turbulence.  Appropriately,  this 
situation  is  referred  to  as  “buoyant  damping.”  The  value  of  can  be  positive  under  the 
ice  if  freezing  is  occurring  at  the  surface.  If  melting  is  occurring  wj’  is  negative 


regardless  of  the  value.  In  the  ice-free  system  is  positive  when  the  net  heat  flux 
is  positive  (out  of  the  system).  When  the  net  heat  flux  is  negative,  or  into  the  ice-free 
system  (summer  conditions),  is  negative  as  well. 
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The  net  effect  of  dissipation  is  absorbed  by  the  coefficients,  c^  and  Cj .  For  this 
work  the  coefficients  Cj  and  Cj  are  considered  tunable  constants.  Kraus  and  Turner 
(1967)  assumed  q  and  Cj  to  be  constant.  First  Gill  and  Turner  (1976),  and  later 
Garwood  (1977)  showed  that  both  q  and  Cj  vary  with  mixed  layer  stability.  However, 
Cj  and  Cj  may  be  considered  constant  for  slowly  changing  surface  forcing  conditions, 

which  are  assumed  in  this  study.  The  parameter  Ha  is  the  "thermobaric  depth" 
(Garwood,  1991)  for  which  increased  pressure  doubles  the  thermal  expansion  relative  to 
the  surface  value,  ,  and  will  be  discussed  further  in  the  following  paragraphs. 

The  buoyancy  flux  is  a  function  of  both  the  potential  temperature  flux  and  the 
salinity  flux.  The  vertical  integral  of  the  buoyancy  flux, 

0  _  0  _  _ 

J b'w\z)dz  =  I  {b'Wg'  +  b'Wg'^dz 

—h  —h 


shows  this  dual  dependence.  Breaking  buoyancy  flux  into  these  two  components  and 
solving  demonstrates  how  each  of  the  fluxes  is  dependent  on  vertical  entrainment.  The 
two  components  are. 
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The  vertical  fluxes  are  related  to  entrainment  by. 


(1)  0'w'-k=-Aew^ 

(2)  'S^'-h=-ASw^ 


where  Ad  and  are  the  jumps  in  potential  temperature  and  salinity  at  the  base  of  the 
mixed  layer.  The  potential  temperature  (0)  and  the  in  situ  temperature  (T)  are  assumed 
to  be  equal  in  this  study,  with  negligible  error.  Taking  the  vertical  integral  of  each 
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component  yields  ,  the  vertical  integral  of  the  buoyancy  flux,  and  its  dependence  on 
the  vertical  entrainment  velocity: 


=  -  PoSh 


S-S, 


(^ogh*  1  +  — 


P.gh*^S 

2 


The  temperature  and  salinity  fluxes  follow  directly  from  the  temperature  and  salinity 
profiles  as  well  as  the  net  heat  flux  into  or  out  of  the  water  column  (Fig.  8).  Temperature 
and  salinity  profdes  in  the  high  latitudes  are  quite  unique,  especially  with  the  presence  of 
ice  (Fig  9).  Assuming  ice  coverage,  the  water  just  below  the  ice  will  be  warmer  than  the 
ice;  if  it  were  not  warmer  it  too  would  be  frozen.  Consequently,  the  water  just  below  the 
ice  is  warmer  than  the  freezing  point  and  usually  continues  to  get  warmer  with  depth. 
Similarly,  since  ice  rejects  salt  when  it  freezes,  the  ice  will  have  a  much  lower  salinity 
than  the  water  beneath  it.  Salinity  usually  continues  to  increase  with  depth.  Recalling 
that  salinity  dominates  over  temperature  in  the  high  latitudes,  these  profiles  create  a 
stable  density  structure. 


When  there  is  ice  coverage  all  of  the  entrainment  works  directly  to  heat  the  ice 
layer  and  does  not  heat  the  mixed  layer  significantly.  When  this  occurs  the  ice  begins  to 
melt  and  so  the  vertical  entrainment  velocity  affects  ice  thickness  (hi)  as  well: 

=  =  03) 

dt  Pi^f  Pi^f 

where  mf  is  the  melting-freezing  term,  quantifies  the  latent  heat  of  fusion,  and  pc^ 
is  the  water  density  times  the  heat  capacity.  When  the  entrainment  heat  flux  exceeds  the 
upward  flux  of  heat  and  radiation  at  the  ice  surface  (Qq),  then  the  ice  will  melt  making 
dk 

<  0  .  Conversely,  when  the  entrainment  is  weaker  than  the  surface  heat  flux  ice  will 


freeze  making  — -  >  0  . 

dt 
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Figure  8.  Representative  temperature  and  salinity  flux  profiles  for  the  high  latitude  regions. 


26 


Atmosphere 


Figure  9.  Representative  salinity  and  temperature  profiles  for  the  high  latitude  regions. 
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In  this  way  vertical  entrainment  simultaneously  affects  both  free  convection  and  the  ice 
thickness.  Therefore,  to  solve  for  entrainment  requires  solving  three  simultaneous 

dh- 

equations,  one  for  free  convection,  wj' ,  one  for  change  in  ice  thickness, - -  ,  and  one 

dt 

for  vertical  entrainment  velocity  itself,  . 

Accuracy  of  the  calculation  of  buoyancy  flux  is  improved  by  considering 
thermobaric  enhancement  of  entrainment.  This  occurs  when  the  thermal  expansion 
coefficient,  a ,  is  no  longer  treated  as  a  constant  but  is  allowed  to  vary  with  both 
temperature  and  depth.  Consequently,  the  equations  for  both  vertical  entrainment 
velocity  and  free  convection  contain  the  term/Za,  the  "thermobaric  depth",  as  well  as  the 
surface  term  for  the  thermal  expansion  coefficient,  .  The  inclusion  of  thermobaric 

depth  can  lead  to  a  magnified  entrainment  rate,  increased  buoyancy,  and  more  efficient 
stirring  of  the  mixed  layer  through  free  convection.  It  is  through  this  process  that  heating 
of  the  ice  can  occur  more  quickly  and  with  greater  intensity.  Prolonged  enhanced 
entrainment  in  a  region  can  lead  to  a  polynya,  where  the  ice  has  been  melted  away 
completely  and  the  entrainment  works  to  warm  the  mixed  layer  to  greater  depth. 

Under  certain  situations  free  convection  can  be  turned  off,  =  0  ,  or  even  be 

made  negative.  As  discussed  previously,  when  this  term  is  negative  it  is  no  longer  called 
free  convection  but  is  instead  called  buoyant  damping.  This  situation  occurs  when  the 
entrainment  shuts  down  or  is  very  weak.  Entrainment,  however,  can  never  be  made 
negative  since  un-stirring,  or  directly  reversing  completed  work  cannot  be  achieved.  If 
free  convection  were  to  shut  down,  forced  convection  could  dominate  and  wind  stirring 
of  the  mixed  layer  could  result.  When  this  happens  the  mixed  layer  could  shallow  to  the 
Obukhov  Length  (L),  and  a  step  would  be  noticed  in  the  density  profile  consistent  with 
summer  shallowing  conditions.  It  is  important  to  note  that  for  the  model  created  in 
association  with  this  work,  circumstances  leading  to  a  negative  sum  of  the  free  and  forced 
convection  will  not  be  processed  accurately. 


D,  MODEL 
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A  one-dimensional  model  eonsisting  of  a  system  of  ordinary  differential  equations 
was  solved  as  an  initial-value  problem  numerieally  to  examine  the  interaetion  between 
the  mixed  layer  and  iee  at  the  surfaee.  The  model  was  designed  this  way  so  that 
partieular  attention  eould  be  paid  to  the  importanee  of  thermobarieity  as  well  as  the 
strength  of  various  foreing  terms  and  initial  eonditions.  Additionally,  this  model  was 
built  with  the  idea  that  it  eould  be  used  direetly  or  as  the  preeursor  for  a  model  that  will 
be  able  to  fit  with  eurrent  iee  predietion  models.  Speeilieally,  the  model  should 
ultimately  work  in  eonjunetion  with  the  Hibler  lee  Model  so  that  the  United  States  Navy 
ean  use  it  operationally  in  systems  like  PIPS  (Fig  .7).  The  end  goal  for  the  future  model 
is  simple;  foreeasts  for  iee  formation,  thiekness,  and  extent  must  all  improve  over  today’s 
model  output.  For  example,  the  iee  thiekness  in  the  Greenland  Sea,  whieh  is  eurrently 
often  poorly  predieted,  may  be  more  aeeurately  represented.  The  model  must  be  robust 
enough  to  generate  iee-free  areas  like  polynyas  or  iee  eovered  regions  given  any  initial 
iee  eondition  and  representative  atmospherie  eonditions.  In  other  words,  given  realistie 
foreing  and  initial  eonditions,  the  model  should  prediet  a  polynya  in  an  area  that  is  iee 
eovered  or  freezing  in  a  polynya.  Matlab’s  C++  was  ehosen  as  the  programming 
language  for  this  model  beeause  of  its  matrix,  differential  equation,  and  graphieal 
eapabilities,  as  well  as  its  ability  to  perform  large  ealeulations  with  minimal  user  effort. 

1,  Model  Equations  Summarized 

a.  Coupled  Mixed  Layer-Ice  System  Equations 

The  system  of  eight  equations  is  derived  by  vertieally  integrating  the 
budgets  for  enthalpy  (heat)(5),  momentum  (2)(3),  salinity  (6)  and  TKE  (8)  between  the 
iee  surfaee  and  the  base  of  the  turbulent  mixed  layer.  Drag  (momentum  flux)  and  heat 
flux  between  the  iee  and  mixed  layer  are  ealeulated  with  bulk  aerodynamie  sealing  laws, 
with  the  temperature  of  the  iee-water  eontaet  assumed  to  be  at  the  freezing  point,  Tfp. 
Entrainment  fluxes  of  heat,  salinity  and  momentum  are  derived  from  integrating  the 
respeetive  eonservation  laws  (2,  3,  5,  and  6)  aeross  the  entrainment  zone  diseontinuity 
(Garwood,  1977). 

Mixed  Layer  Depth; 


29 


dh  I  3 

-r  =  w^={cu,  -c. 
dt 


ag-^  + Pg{mf)({S)-S,) 

p 


■l{agM-PgKS]  (14) 


Mixed-Layer  Temperature; 
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Mixed-Layer  Salinity: 
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Momentum  and  Heat  Flux  Between  Ice  and  Mixed  Layer: 


T,-PCD(ur^)J{ui-u)  +(v/-v)  (22) 


Ty.-PCD(vr^}J(ui—^)  +(v/-v)  (23) 


30 


T'w' 


(24) 


2.  Numerical  Methods 

Matlab  is  a  programming  language  that  can  perform  large  calculations  and 
provide  graphical  representation  of  data  through  the  manipulation  of  matrices  and 
execution  of  user  produced  or  existing  functions.  Its  ability  to  quickly  and  accurately 
perform  complicated  and  or  lengthy  computations  make  it  an  attractive  tool  for  ocean 
modeling. 

a.  Use  of  Stiff  Solver 

ODE23s  is  a  Matlab  function  that  solves  ordinary  differential  equations 
using  a  stiff  solver.  The  use  of  ODE23s  allowed  for  the  tracking  and  graphing  of  many 
different  oceanographic  parameters  as  the  program  stepped  forward  in  time.  Specifically, 
it  kept  track  of  and  depicted  how  both  the  mixed  layer  depth  and  ice  thickness  increased 
or  decreased  given  the  initial  and  forcing  conditions.  The  program  ODE23s  was  chosen 
over  other  Matlab  differential  equation  solvers  because  of  its  ability  to  solve  stiff 
differential  equations.  A  way  of  solving  this  type  of  equation  became  necessary  when 
approaching  or  crossing  into  the  ice-free  regime,  where  intermediate  solutions  change 
rapidly,  but  the  final  solution  is  found  on  a  longer  time  scale  (Shampine,  1994).  Not 
using  a  stiff  solver  would  mean  allowing  for  the  time  step  to  decrease  so  that  every 
intermediate  solution  could  be  calculated.  This  ability  was  not  necessary  since  the  final 
solution  was  the  solution  of  interest.  To  numerically  solve  these  types  of  equations 
Matlab  employs  the  Runge-Kutta  formula  or,  in  the  case  of  ODE23s,  a  modified 
Rosenbrock  formula  that  is  more  efficient  (Math  Works,  2002).  Additionally,  because  of 
the  computational  power  and  time  benefits  of  this  program,  it  was  verified  that  treating 
the  water  column  as  three  slabs  was  a  robust  method  for  the  desired  modeling.  This  was 
proven  graphically  by  showing  that  the  total  of  the  velocities  of  the  ice  layer  slab  and 
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mixed  layer  slab,  individually,  were  equal  in  value  to  the  veloeity  of  a  combined  ice- 
mixed  layer  slab. 
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IV.  SENSITIVITY  STUDIES 


The  following  sensitivity  studies  simulated  Weddell  Sea  winter  conditions. 
Akitomo  (1999)  first  applied  Garwood  et  al.  (1994)  thermobaric  instabilities  to  the 
Weddell  Sea,  but  without  the  ice-mixed  layer  feedback  model  used  here.  This  was  done 
because  there  is  greater  interaction  between  temperature  and  salinity  in  this  region  and  at 
this  time.  For  each  simulation  one  of  the  forcing  conditions,  net  heat  flux  (Qq)  or  wind 

speed,  was  held  constant  while  the  other  was  allowed  to  vary.  Similarly,  the  initial 
conditions  were  varied  for  the  simulations.  Two  different  depths  of  the  mixed  layer  were 
chosen  to  examine  the  effects  of  the  ratio  between  mixed  layer  depth  and  thermobaric 
depth.  These  depths  were  100  m  for  a  moderately  deep  mixed  layer  and  1000  m  for  a 
deep  mixed  layer.  The  effects  of  changing  the  ratio  of  thermal  stratification  to  salinity 
stratification,  a^gAT  to  /^ggAS ,  on  ice  melting  and  freezing  at  the  surface  were 
examined. 

It  should  be  noted  that  the  program  allows  for  the  ratio  of  CC^gAT  to  P^gAS  to 
change  with  changing  temperature  and  salinity  values.  If  these  values  were  instead 
prescribed,  accuracy  would  be  lost  since  it  would  effectively  not  allow  for  any  storage  of 
heat  or  salt  in  the  mixed  layer  when  ice  is  present  at  the  surface.  Furthermore,  any  warm 
salty  water  entrained  would  not  change  the  mixed  layer  temperature  and  salinity  values. 
Equations  from  the  previous  chapter  demonstrate  how  important  the  temperature  and 
salinity  fluxes  within  the  mixed  layer  are  to  the  free  convection  and  total  turbulent  kinetic 
energy  of  the  system.  Not  allowing  for  these  fluxes  is  the  equivalent  of  neglecting 
physical  processes  essential  to  the  coupled  ice-mixed  layer  system.  Instead,  by  allowing 
for  this  change  in  the  ratio,  the  vertical  integral  of  the  buoyancy  flux  can  be  accurately 
calculated. 

In  all  simulations  the  initial  ice  thickness  is  10  cm;  the  initial  temperature  is  the 
temperature  of  the  freezing  point  at  the  surface,  -1.66  Celsius;  the  initial  salinity  value  is 
34.52  parts  per  thousand;  and  the  thermobaric  depth  is  93.6  m.  Since  both  100  m  and 
1 000  m  mixed  layers  are  deeper  than  the  thermobaric  depth,  thermobaricity  is  expected  to 
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enhance  the  entrainment  of  warmer  and  saltier  water  from  below.  Each  case  is  run  for 
ten  days,  or  8.64x10^  seconds. 

A,  MIXED  LAYER  DEPTH  EQUAL  TO  100  METERS 

For  the  first  simulation  at  h=100  m,  the  ratio  of  cx^gAT  to  ^^gAS  equaled  0.45  for 

the  surface  temperature  and  salinity  values.  The  total  heat  flux  was  a  constant  0.001 
cal/cm  s  upward  and  out  of  the  coupled  system.  The  wind  speed  was  allowed  to  vary 
from  400  cm/s  to  1600  cm/s  in  400  cm/s  increments.  The  results  show  that  for  each  wind 
speed  positive  entrainment  deepens  the  mixed  layer  (Fig.  10).  The  weakest  wind,  400 
cm/s,  produced  the  least  amount  of  mixed  layer  deepening.  The  resultant  entrainment 
was  not  strong  enough  to  melt  ice  at  the  surface  and  instead  fostered  ice  growth. 

The  mixed  layer  salinity  increases  for  the  400  cm/s  wind  as  freezing  occurs  and 
salt  is  rejected  from  the  ice.  The  corresponding  mixed  layer  temperature  warms  until  it 
comes  into  equilibrium.  Once  the  temperature  equilibrium  is  reached,  the  rate  of  salinity 
increase,  ice  growth,  and  mixed  layer  deepening  all  become  constant. 

The  higher  wind  values  of  800  cm/s,  1200  cm/s,  and  1600  cm/s  each  produced 
strong  enough  entrainment  to  melt  the  ice  and  deepen  the  mixed  layer,  with  the  1600 
cm/s  wind  having  the  strongest  melting  rate.  Similarly,  the  1600  cm/s  case  had  the 
largest  decrease  in  salinity  associated  with  the  melting  of  the  ice  and  injection  of  fresh 
water  into  the  mixed  layer.  Correspondingly,  the  mixed  layer  temperature  warmed  the 
most  significantly  immediately  before  eventually  approaching  equilibrium.  It  should  be 
noted  that  the  1600  cm/s  plots  do  not  extend  to  the  full  ten  days  like  the  other  wind 
speeds.  This  is  due  to  the  strong  entrainment  that  occurs  and  the  ill-defined  matrix  that 
results  when  trying  to  continue  to  solve  for  the  entrainment  and  terms 
simultaneously. 
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Figure  10.  Simulation  #1:  Sensitivity  study  for  varying  winds  with  h  =100  m,  Qo  =0.0010 
eal/cm^s,  and  ratio  of  Gf^gAT  to  /3QgAS  equal  to  0.45.  Wind  speed  ehanges  from  400 
ern/s  to  1600  em/s  in  400  cm/s  intervals. 
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Winds  of  1200  cm/s  and  1600  cm/s  are  capable  of  melting  away  the  iee 
eompletely;  however,  onee  the  iee  is  melted  away  the  model  shows  negative  iee  thiekness 
at  approximately  2.3  days  (2.0x10^  see)  and  6.9  days  (6. 0x10^ see)  respeetively.  This  is 
an  unrealistie  depietion  of  the  true  physieal  proeess;  there  is  no  sueh  thing  as  negative  iee. 
Similarly,  an  unrealistie  deerease  in  salinity  values  is  aehieved  for  these  wind  speeds  at 
the  same  points  in  time.  The  salinity  values  are  unrealistie  beeause  onee  the  iee  has 
melted  away  there  is  no  other  souree  of  fresh  water  to  eontinue  to  deerease  the  salinity  of 
the  mixed  layer.  Future  versions  of  the  model  need  to  be  able  to  melt  away  the  iee, 
remain  iee-free  without  predieting  negative  iee,  and  then  be  able  to  make  iee,  beginning 
with  no  ice,  when  eonditions  warrant  iee  growth.  For  this  situation  the  iee  should  melt 
away  due  to  strong  entrainment  and  then  the  system  should  remain  iee  free  due  to  a 
sustained  high  value  of  entrainment.  The  high  entrainment  rate  should  eontinue 
deepening  even  without  iee  at  the  surfaee,  although  the  true  deepening  rate  is  probably 
different  than  the  one  represented  here  with  negative  iee  thiekness. 

Examining  the  plot  of  mixed  layer  temperature  with  time  reveals  that  there  is 
more  heat  storage  in  the  mixed  layer  beneath  freezing  eonditions  (wind  speed  equals  400 
em/s)  at  the  surfaee  than  there  is  with  melting  eonditions  (800-1600  em/s  wind  speeds)  at 
the  surfaee.  That  is,  the  mixed  layer  equilibrium  temperature  is  warmer  beneath  iee 
growing  eonditions  than  it  is  beneath  iee  melting  eonditions.  This  makes  sense  when  the 
iee  at  the  surfaee  is  thought  of  as  a  eap  that  does  not  allow  the  entrained  heat  to  eseape 
out  of  the  system.  Sinee  it  is  not  strong  enough  to  melt  the  iee  and  eannot  eseape  to  the 
atmosphere,  this  heat  is  stored  in  the  mixed  layer.  Also  shown  through  this  sensitivity 
study  is  that  for  a  given  net  heat  flux,  wind  ean  eause  entrainment  to  be  strong  enough  to 
melt  iee  even  in  freezing  atmospherie  eonditions.  Furthermore,  strong  entrainment  ean 
eause  and  sustain  polynya-like  eonditions. 

For  the  seeond  simulation  at  h=100  m,  the  ratio  of  (X^gAT  to  P^gAS  equaled 
0.036  for  the  surfaee  temperature  and  salinity  values.  The  total  heat  flux  was  a  eonstant 
0.001  eal/em  s  upward  and  out  of  the  eoupled  system.  The  wind  speed  was  allowed  to 
vary  from  400  em/s  to  1600  em/s  in  400  em/s  inerements.  The  results  of  this  wind 
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sensitivity  study  show  how  important  the  ratio  of  cx^gAT  to  J^^gAS  is  to  entrainment  and 
ice  formation  or  melting. 

Like  the  first  simulation,  the  mixed  layer  depth  increases  for  all  wind  speeds  due 
to  positive  entrainment  (Fig.l  1).  The  1600  cm/s  wind  caused  stronger  entrainment  than 
the  400  cm/s  wind  as  well.  However,  for  this  case  the  entrainment  is  never  strong  enough 
to  melt  the  ice.  Entrainment  strength  is  still  dependent  on  the  wind  speed,  with  the  higher 
wind  speeds  creating  greater  entrainment.  This  can  be  seen  in  the  freezing  rate  where  the 
lower  wind  speeds  with  weaker  entrainment  have  higher  freezing  rates.  In  this  case  the 
400  cm/s  wind  freezes  at  a  rate  of  approximately  1.1  cm/day  (11  cm  in  ten  days)  as 
compared  to  the  1600  cm/s  wind  that  freezes  at  a  rate  of  approximately  0.4  cm/day  (4  cm 
in  ten  days). 

Examining  the  ice  thickness  time  series  plot  shows  that  the  1600  cm/s  wind  slows 
the  freezing  rate.  Given  a  few  more  days  with  this  constant  wind  speed  it  is  conceivable 
that  the  ice  formation  would  halt  and  that  melting  might  even  begin.  Salinity  predictably 
increases  as  freezing  occurs  at  all  wind  speeds.  Here,  however,  the  salinity  increases  the 
most  for  the  1600  cm/s  wind  as  compared  to  lesser  winds.  This  is  in  contrast  to  the  result 
from  the  first  simulation  where  the  greatest  freezing  rate  had  the  greatest  increase  in 
mixed  layer  salinity.  Comparing  the  400  cm/s  wind  to  the  1600  cm/s  wind  and  the 
resulting  ice  formation  of  both  it  is  seen  that  the  400  cm/s  wind  produces  a  weaker 
entrainment,  greater  ice  thickness,  and  a  less  saline  mixed  layer  as  compared  to  the  1600 
cm/s  wind  that  produces  greater  entrainment,  thinner  ice,  and  a  more  saline  mixed  layer. 
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Figure  1 1 .  Simulation  #2;  Sensitivity  study  for  varying  winds  with  h  =100  m,  Qo  =0.0010 
cal/em^s,  and  ratio  of  CC^g/ST  to  ^^gAS  equal  to  0.036.  Wind  speed  changes  from  400 

cm/s  to  1600  cm/s  in  400  cm/s  intervals. 
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The  400  cm/s  wind  is  the  slowest  to  bring  the  system  into  thermal  equilibrium, 
but  appears  to  achieve  a  higher  equilibrium  temperature  than  some  of  the  greater  wind 
speeds.  This  can  be  seen  by  looking  at  the  mixed  layer  temperature  plot  at  day  ten  where 
the  400  cm/s  wind  has  already  surpassed  the  mixed  layer  temperature  value  for  the  800 
cm/s  wind  and  appears  close  to  exceeding  the  mixed  layer  temperature  of  the  1200  cm/s 
wind.  Like  in  the  first  simulation,  where  ice  is  present  and  freezing  is  occurring,  the 
mixed  layer  has  a  higher  equilibrium  temperature  leading  to  more  storage  of  heat  from 
the  effect  of  ice  capping  the  system.  These  results  indicate  that  strong  entrainment  with 
ice  presence  will  bring  the  system  into  thermal  equilibrium  faster,  but  that  entrainment  is 
not  the  only  factor  for  establishing  what  that  point  of  thermal  equilibrium  will  be. 

For  the  third  simulation  at  h=100  m,  the  ratio  of  a^gAT  to  /^^gAS  again  equaled 

0.45  for  the  surface  temperature  and  salinity  values.  The  wind  speed  was  a  constant 
value  of  700  cm/s.  The  heat  flux  was  allowed  to  vary  from  0.0005  caFcm  s  to  0.0025 
cal/cm  s  in  0.0005  cal/cm  s  increments.  The  results  show  that  for  all  values, 

sufficient  entrainment  was  produced  to  deepen  the  mixed  layer  (Fig.  12).  Similar  to  the 
second  simulation,  the  stronger  the  heat  flux  the  greater  entrainment  it  caused  and  the 
deeper  the  mixed  layer  grew.  A  low  value  equal  to  0.0005  cal/cm  s  causes  deepening 

to  approximately  103  m,  while  a  value  equal  to  0.0025  cal/cm  s  causes  deepening  to 
1 15  m,  for  a  difference  of  12  meters  of  deepening  between  the  two  fluxes. 

For  all  cases  of  gg  the  ice  actually  grows  at  first  but  then  begins  to  melt.  The 
time  that  the  transition  from  freezing  to  melting  occurs  differs  for  each  value; 
however,  the  greater  the  Qq  value  the  longer  it  takes  for  this  transition  to  occur.  For 
example,  for  a  heat  flux,  gg  equal  to  0.0005  cal/cm  s  it  takes  approximately  1.2  days 
(1.0x10  sec)  for  melting  to  begin.  A  heat  flux,  ,  equal  to  0.0025  cal/cm  s  will  take  4.0 
days  (3.5x10^  sec)  for  this  change  from  freezing  to  melting  to  occur. 
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Figure  12.  Simulation  #3;  Sensitivity  study  for  heat  fluxes  with  h  =100  m,  wind  speed  =700 
cm/s,  and  ratio  of  (Xf^gAT  to  /^QgAS  equal  to  0.45.  Heat  flux  changes  from  0.0005 
cal/cm^s  to  0.0025  cal/cm^s  in  0.0005  cal/cm^s  intervals. 
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Mixed  layer  salinity  inereases  when  ice  is  growing  and  decreases  when  ice  is 
melting.  This  is  consistent  with  the  previous  two  simulations  and  the  physics  of  salt 
rejection  during  ice  freezing.  When  equal  to  0.0005  cal/cm  s  begins  melting  ice,  the 
salinity  predictably  starts  to  decrease,  since  the  melting  puts  fresh  water  into  the  mixed 
layer.  However,  at  larger  values  of  gg  when  melting  occurs  the  salinity  does  not 

decrease  as  it  does  with  the  equal  to  0.0005  cal/cm  s;  instead,  the  rate  of  salinity 
increase  becomes  constant.  The  explanation  for  this  is  that  entrainment  is  lifting  more 
salinity  into  the  mixed  layer  from  below  than  it  is  taking  out  of  the  mixed  layer  through 
melting  and  fresh  water  addition  to  the  mixed  layer. 

High  heat  flux  values,  like  0.0025  cal/cm  s,  that  lead  to  relatively  strong  freezing 
had  the  highest  mixed  layer  equilibrium  temperature,  and  also  approached  equilibrium 
faster  than  lower  heat  flux  values.  These  results  can  be  observed  by  comparing  equal 
to  0.0005  cal/cm^s  with  equal  to  0.0020  cal/cm^s.  For  0.0005  cal/cm^s  the  mixed 
layer  temperature  came  to  equilibrium  at  approximately  -1.64°C  after  3.5  days 
(3.0x10  sec),  while  the  0.0020  cal/cm  s  heat  flux  came  to  equilibrium  at  approximately  - 
1 .597C  after  8.7  days  ( 7.5x10^  sec).  For  these  measurements  the  rate  at  which  the  mixed 
layer  temperature  comes  to  equilibrium  for  the  lower  heat  flux  is  .0057C/day  and  the  rate 
at  which  the  mixed  layer  temperature  comes  to  equilibrium  for  the  higher  valued  heat 
flux  is  .0069C/day.  These  results  suggest  that  there  may  be  a  relationship  between  ice 
thickness  and  heat  storage.  This  relationship  would  come  into  practice  when  strong 
entrainment  does  little  melting  but  mixed  layer  depth  continues  to  increase  (deepen). 
Here,  where  the  ice  is  not  melted  and  heat  cannot  escape  to  the  atmosphere,  the  heat 
becomes  stored  in  the  mixed  layer  and  causes  warming  of  the  equilibrium  temperature. 

Comparing  this  simulation  to  the  first  simulation,  where  the  heat  flux  was  kept 
constant,  produces  some  interesting  results.  First,  the  lowest  Qg  value  produces  the 

coolest  equilibrium  temperature,  whereas  the  lowest  wind  speed  appeared  to  have  the 
warmest  equilibrium  temperature.  This  indicates  that  more  heat  can  be  stored  in  the 
mixed  layer  of  the  first  simulation  as  compared  to  the  mixed  layer  of  the  third  simulation 
for  minimum  forcing  conditions.  Second,  a  400  cm/s  wind  with  constant  Qg  produced 
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freezing  of  ice,  but  the  lowest  value,  0.0005  cal/cm  s,  with  constant  wind  gives 

melting.  Both  of  these  observations  can  be  explained  by  concluding  that  the  constant 
values  used  in  simulation  are  weak  enough  to  allow  other  factors  to  influence  the 
entrainment  and  thus  change  the  characteristics  of  the  coupled  system.  The  usefulness  of 
these  constant  values  for  heat  flux  and  wind  speed  is  further  substantiated  because  in 
simulation  1  where  gg  is  constant  and  winds  are  weak,  freezing  is  forecast.  Similarly,  in 

simulation  3  where  wind  speed  is  constant  and  gg  is  weak,  melting  is  forecast. 

For  the  fourth  simulation  at  h=100  m,  the  ratio  of  to  /^^gAS  equaled  0.036 

for  the  surface  temperature  and  salinity  values.  The  wind  speed  was  a  constant  value  of 
700  cm/s.  The  heat  flux  was  allowed  to  vary  from  0.0005  cal/cm  s  to  0.0025  cal/cm  s  in 
0.0005  cal/cm  s  increments.  The  results  demonstrate  that  positive  entrainment  causes 
mixed  layer  deepening  for  all  heat  flux  values  (Fig.  13).  The  ice  thickness  growth  rate 
increases  as  increases  from  0.0005  cal/cm  s  to  0.0025  cal/cm  s  with  no  transition  to 
melting.  Salinity  predictably  increases  for  all  heat  flux  values,  with  the  strongest  rate  of 
salinity  increase  corresponding  to  equal  to  0.0025  cal/cm  s.  Like  simulation  3,  the 

greatest  change  in  ice  thickness  is  associated  with  the  warmest  equilibrium  temperature 
and  the  largest  amount  of  heat  storage  within  the  mixed  layer. 
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Figure  13.  Simulation  #4:  Sensitivity  study  for  heat  fluxes  with  h  =100  m,  wind  speed  =700 
cm/s,  and  ratio  of  a^gAT  to  /3f^gAS  equal  to  0.036.  Fleat  flux  changes  from  0.0005 
cal/cm^s  to  0.0025  cal/cm^s  in  0.0005  cal/cm^s  intervals. 
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When  comparing  this  simulation  with  the  previous  one  the  effects  of  the  different 
ratios  of  cx^gAT  to  ^ogAS  can  be  examined.  It  appears  as  though  the  smaller  ratio, 
0.036,  increases  the  significance  of  while  simultaneously  weakening  the  role  of 
entrainment.  A  weakening  of  entrainment  would  explain  why  the  mixed  layer  does  not 
deepen  as  much  in  simulation  4  as  it  does  in  simulation  3.  The  strongest  value  for  in 
the  third  simulation  deepens  the  mixed  layer  to  approximately  115  m.  The  strongest 
value  for  in  the  fourth  simulation  only  deepens  the  mixed  layer  to  103.5  m.  Similarly, 

the  strongest  value  of  Qq  in  the  third  simulation  produces  only  approximately  1 1.5  cm  of 
ice,  while  the  strongest  value  of  in  the  fourth  simulation  produces  approximately  37 

cm  of  ice.  The  fact  that  the  mixed  layer  temperature  remains  cooler  for  each  heat  flux 
value  in  simulation  4  as  compared  to  the  same  heat  flux  values  in  simulation  3  lends 
credence  to  the  possibility  that  entrainment  has  weakened,  since  less  heat  is  stored  in  the 
mixed  layer  of  simulation  4.  Further  evidence  of  weakened  entrainment  is  found  in  the 
mixed  layer  salinity  versus  time  plot  where  for  simulation  4  the  salinity  values  for  each 
heat  flux  value  reach  higher  values  than  they  do  in  the  third  simulation.  This  corresponds 
to  the  net  effect  of  increased  salt  rejection  and  no  addition  of  fresh  water  from  melting 
caused  by  entrained  heat. 

B.  MIXED  LAYER  DEPTH  EQUAL  TO  1000  METERS 

For  the  first  simulation  at  h=1000  m,  the  ratio  of  a^gAT  to  P^gAS  equaled  0.45 
for  the  surface  temperature  and  salinity  values.  The  total  heat  flux  was  a  constant  0.001 
cal/cm  s  up  and  out  of  the  coupled  system.  The  wind  speed  was  allowed  to  vary  from 
400  cm/s  to  1600  cm/s  in  400  cm/s  increments.  The  results  show  that  positive 
entrainment  deepens  the  mixed  layer  for  all  of  the  wind  speeds  (Fig.  14).  The  400  cm/s 
wind  causes  sustained  deepening  by  creating  relatively  weak  entrainment  that  is  not 
strong  enough  to  melt  the  ice.  The  three  other  larger  values  for  wind  cause  freezing 
initially  but  then  shift  to  melting  conditions  when  entrainment  overpowers  the  net  heat 
flux.  This  initial  freezing  period  lasts  approximately  2.3  days  ( 2.0x10^  sec).  After  this 
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period,  the  mixed  layer  deepens  twiee  as  fast  in  the  400  cm/s  case  as  it  does  in  the  800 
cm/s  case. 

This  rate  can  be  determined  by  picking  values  off  of  the  mixed  layer  depth  versus 
time  plot  for  both  the  400  cm/s  line  and  800  cm/s  line  corresponding  to  times  of 
2.0x10^  seconds  and  4.0x10^  seconds  (4.6  days).  For  the  400  cm/s  wind  at  time  equal  to 
2.0x10^  seconds  the  mixed  layer  depth  equals  1013  m,  and  at  time  equal  to 
4.0x10^  seconds  the  mixed  layer  depth  equals  1021  m.  For  the  800  cm/s  wind  at  time 
equal  to  2.0x10^  seconds  the  mixed  layer  depth  equals  1010  m,  and  at  time  equal  to 
4.0x10^  seconds  the  mixed  layer  depth  equals  1014  m.  Solving  for  the  change  in  depth 
per  day  for  each  case  yields  a  change  of  3.47  m/day  for  the  400  cm/s  wind  and  a  change 
of  1.74  m/day  for  the  800  cm/s  wind. 

The  ice  thickness  versus  time  plot  changes  according  to  the  strength  of 
entrainment  where  greater  entrainment  causes  melting  to  occur  sooner.  Predictably, 
salinity  increases  with  increasing  freezing  and  decreases  for  strong  melting,  as  is  the  case 
with  the  1600  cm/s  wind.  The  400  cm/s  wind  has  the  highest  equilibrium  temperature  of 
the  four  wind  speeds.  Conversely,  the  1600  cm/s  wind  has  the  lowest  equilibrium 
temperature  but  also  achieves  thermal  equilibrium  more  quickly  than  the  lesser  wind 
speeds.  This  result  is  due  to  the  fact  that  the  high  winds  produce  strong  entrainment  that 
is  enhanced  by  thermobaricity  to  melt  the  ice  without  significantly  warming  the  mixed 
layer.  Like  thermobaricity,  the  ratio  of  mixed  layer  depth  to  thermobaric  depth  likely 
plays  a  role  in  the  sustained  deepening  of  the  mixed  layer  for  the  400  cm/s  wind.  The 
ration  of  the  mixed  layer  to  the  thermobaric  depth  was  calculated  using  depth  of  the 
mixed  layer,  equal  to  1000  m,  divided  by  thermobaric  depth,  which  is  equal  to  93.6  m  for 
all  of  the  simulations. 
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Figure  14.  Simulation  #5:  Sensitivity  study  for  varying  winds  with  h  =1000  m,  Qo  =0.0010 
cal/cm^s,  and  ratio  of  Gf^gAT  to  /3QgAS  equal  to  0.45.  Wind  speed  ehanges  from  400 
cm/s  to  1600  cm/s  in  400  cm/s  intervals. 
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For  the  second  simulation  at  h=1000  m,  the  ratio  of  (X^gAT  to  P^gAS  equaled 

0.036  cal/cm  s  for  the  surface  temperature  and  salinity  values.  The  total  heat  flux  was  a 
constant  0.001  cal/cm^s  up  and  out  of  the  coupled  system.  The  wind  speed  was  allowed 
to  vary  from  400  cm/s  to  1600  cm/s  in  400  cm/s  increments.  The  results  show  that 
positive  entrainment  causes  deepening  for  all  cases,  where  the  stronger  the  wind  was  the 
greater  the  mixed  layer  deepened  (Fig.  15).  For  this  simulation  entrainment  never 
became  strong  enough  for  melting  to  occur  as  evidenced  by  the  constant  rates  of 
deepening  for  all  wind  speeds.  The  freezing  rate  was  strongest  for  the  lightest  wind  but 
all  of  the  wind  speeds  had  very  similar  freezing  rates.  Not  surprisingly  increase  in 
salinity  for  all  wind  speeds  was  observed. 

Looking  closely  at  the  temperature  scale  on  the  y-axis  of  the  mixed  layer 
temperature  versus  time  plot  shows  that  the  mixed  layer  temperature  began  near 
equilibrium  and  remained  relatively  constant  for  this  time  period.  Since  there  is  very 
little  heat  storage  in  the  mixed  layer  with  ice  present  at  the  surface  it  seems  possible  that 
entrainment  is  either  very  weak  or  has  not  yet  become  significant  as  compared  to  the  net 
heat  flux. 

As  compared  to  the  previous  simulation  this  simulation’s  mixed  layer  does  not 
deepen  as  much  for  any  of  the  wind  speeds.  Similarly,  the  mixed  layer  remains  cooler 
and  very  little  heat  storage  takes  place  for  all  the  wind  speeds  in  this  simulation  as 
compared  to  the  fifth  simulation.  On  average  the  mixed  layer  salinity  increases  more  for 
the  sixth  simulation  than  it  does  for  the  fifth  simulation.  However,  for  the  400  cm/s  case 
the  salinity  only  reaches  34.525  ppt  as  compared  to  the  slightly  more  saline  value  of 
34.526  ppt  from  simulation  5. 
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Figure  15.  Simulation  #6;  Sensitivity  study  for  varying  winds  with  h  =1000  m,  Qo  =0.0010 
cal/em^s,  and  ratio  of  a^g/ST  to  ^^gAS  equal  to  0.036.  Wind  speed  ehanges  from  400 

cm/s  to  1600  cm/s  in  400  cm/s  intervals. 
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For  the  third  simulation  at  h=1000  m,  the  ratio  of  (XffgAT  to  equaled  0.45  for  the 

surfaee  temperature  and  salinity  values.  The  wind  speed  was  a  eonstant  value  of  700 
em/s.  The  heat  flux  was  allowed  to  vary  from  0.0005  eahem^s  to  0.0025  eal/em^s  in 
0.0005  eal/em  s  inerements.  The  results  show  that  positive  entrainment  eauses  deepening 
of  the  mixed  layer  for  all  values  of  Qq  (Fig.  16).  There  is  a  greater  entrainment  for  the 
higher  values  of  heat  flux  demonstrated  by  the  inereased  deepening  of  the  mixed  layer  for 
the  higher  values  of  heat  flux.  For  example,  for  equal  to  0.0025  eal/em  s  the  mixed 

layer  depth  inereases  to  approximately  1055  m,  whereas  the  mixed  layer  depth  for  Qq 
equal  to  0.0005  eal/em  s  only  inereases  to  approximately  1010  m.  For  all  values  of 
entrainment  causes  a  shift  from  freezing  to  melting  conditions.  The  heat  flux  equal  to 
0.0005  cal/cm^s  transitions  first  from  freezing  to  melting  at  day  4.6  (4.0x10^  sec).  Each 
subsequent  heat  flux  transitions  later  than  its  predecessor  until  gg  equal  to  0.0025 

cal/cm^s  transitions  at  day  5.2  ( 4.5x10^  sec).  This  transition  occurs  more  quickly  for 
lower  heat  flux  values  because  even  though  the  entrainment  produced  by  this  low  heat 
flux  value  is  less  than  the  entrainment  produced  by  higher  heat  flux  values,  it  is  more 
easily  overcome  by  the  entrainment. 

Salinity  follows  predictably  from  the  mixed  layer  depth  plot.  Again,  even  though 
melting  occurs  the  salinity  continues  to  increase  due  to  entrainment  of  higher  salinity 
water  from  below.  Significant  heat  storage  in  the  mixed  layer  can  be  seen  where  the 
mixed  layer  equilibrium  temperature  is  warmest  for  the  highest  value.  As  before,  the 

highest  equilibrium  temperature  appears  linked  to  the  highest  ice  thickness  so  that 
entrainment  below  ice  brings  heat  up  into  the  mixed  layer  and  holds  it  there  while  the 
entrained  heat  works  to  melt  the  ice. 
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Figure  16.  Simulation  #7:  Sensitivity  study  for  heat  fluxes  with  h  =1000  m,  wind  speed  =700 
cm/s,  and  ratio  of  to  equal  to  0.45.  Heat  flux  changes  from  0.0005 

cal/cm^s  to  0.0025  cal/cm^s  in  0.0005  cal/cm^s  intervals. 
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Comparing  this  simulation  with  simulation  5  so  that  the  effects  of  keeping  gg 

constant  versus  keeping  wind  speed  constant  produce  interesting  results.  For  this 
simulation  the  mixed  layer  depth  can  increase  by  as  much  as  25  m  over  the  maximum 
mixed  layer  depth  of  simulation  5.  However,  in  simulation  5  the  maximum  depth  was 
produced  by  the  400  cm/s  wind,  whereas  in  this  simulation  the  deepest  mixed  layer  was 
produced  by  the  0.0025  caFcm  s  heat  flux.  This  demonstrates  that  the  constant  values  for 
Qq  and  wind  speed  employed  are  weak  enough  to  allow  other  factors  to  influence  the 
entrainment  and  thus  change  the  characteristics  of  the  coupled  system. 

Comparing  this  simulation  with  the  third  simulation  shows  the  effects  of  differing 
initial  mixed  layers.  The  effects  include;  increased  change  of  mixed  layer  depth  from 
approximately  15  m  for  the  greatest  heat  flux  in  simulation  3,  to  55  m  for  the  strongest 
heat  flux  in  simulation  7;  increased  maximum  salinity  value  for  the  greatest  heat  flux  in 
simulation  3  (34.547  ppt),  as  compared  to  the  maximum  salinity  value  for  the  greatest 
heat  flux  in  simulation  7  (34.530  ppt);  and  mixed  layer  equilibrium  temperature  while 
relatively  constant  for  all  of  the  heat  flux  values  between  the  two  simulations  does 
increase  faster  in  simulation  3  than  in  simulation  7,  where  equilibrium  is  reached  in 
simulation  3  after  approximately  3.5  days  (3.0x10^  sec)  as  compared  to  8.7  days 
(7.5x10^  sec)  for  simulation  7. 

For  the  fourth  simulation  at  h=1000  m,  the  ratio  of  (X^gAT  to  P^gAS  equaled 

0.036  for  the  surface  temperature  and  salinity  values.  The  wind  speed  was  a  constant 
value  of  700  cm/s.  The  heat  flux  was  allowed  to  vary  from  0.0005  caFcm  s  to  0.0025 
cal/cm^s  in  0.0005  cal/cm^s  increments.  The  results  show  that  positive  entrainment 
creates  mixed  layer  deepening  for  all  heat  flux  values  (Fig.  17).  Entrainment  is  weak 
since  the  mixed  layer  does  not  increase  much  with  time  (maximum  of  3.5  m),  freezing 
occurs  for  all  heat  flux  cases  with  the  maximum  ice  thickness  equal  to  39  cm  for 

equal  to  0.0025  cal/cm  s,  and  very  little  heat  can  be  stored  in  the  mixed  layer  since  the 
mixed  layer  temperature  varies  little  even  in  the  most  extreme  heat  flux  case  (0.003C 
over  10  days). 
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Comparing  this  simulation  with  simulation  4  (h=100  m)  to  examine  the  effeets  of 
initial  mixed  layer  depth  produees  the  following:  3.5  meters  of  depth  inerease  for  both 
simulations  where  equals  to  0.0025  eal/em  s;  simulation  8  produees  approximately  2 

em  more  iee  for  all  heat  flux  eases  than  simulation  4;  the  mixed  layer  temperature  warms 
faster  and  to  larger  equilibrium  values  in  simulation  4  than  it  does  in  simulation  8;  and 
the  mixed  layer  salinity  inereases  faster  and  to  higher  values  in  simulation  4  than  it  does 
in  simulation  8.  These  results  ean  be  verified  by  looking  at  the  slope  of  the  respeetive 
lines  or  by  taking  sample  points  off  of  the  lines  and  eomparing  the  change  in  a  quality 
over  time  numerically  as  done  on  several  occasions  above. 

Comparing  this  simulation  with  simulation  6  to  examine  the  effects  of  gg  versus 

wind  forcing  on  the  deep  mixed  layer  indicates  that  the  effect  of  strong  heat  flux  values 
out  of  the  system  is  more  significant  to  freezing  than  low  wind  speeds.  This  can  be 
quantified  by  calculating  the  ice  growth  from  simulation  6  for  low  wind  (400  cm/s)  to  be 
equal  to  12  cm.  Similarly,  calculating  the  ice  growth  from  simulation  8  for  high  heat  flux 
(0.0025  cal/cm  s)  yields  30  cm  of  ice.  From  these  results  it  can  be  concluded  that  for 
these  values  of  heat  flux  and  wind  speed,  large  heat  flux  tends  to  produce  stronger 
freezing  conditions  than  low  wind  forcing. 
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Figure  17.  Simulation  #8;  Sensitivity  study  for  heat  fluxes  with  h  =1000  m,  wind  speed  =700 
cm/s,  and  ratio  of  a^gAT  to  /^^gAS  equal  to  0.036.  Fleat  flux  changes  from  0.0005 
cal/cm^s  to  0.0025  cal/cm^s  in  0.0005  cal/cm^s  intervals. 
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Comparing  this  simulation  with  simulation  7  to  examine  the  effeets  of  the  ratio  of 
a^gAT  to  /^^gAS  shows  that  the  higher  ratio,  like  that  of  simulation  7,  produees  melting 

of  the  ice  and  greater  deepening  of  the  mixed  layer  (up  to  55  m  gained  as  compared  to  a 
maximum  of  32  m  gained  in  simulation  8).  Simulation  7  also  can  have  greater  heat 
storage  within  the  mixed  layer  as  a  result  of  the  increased  warming  of  the  mixed  layer. 
The  mixed  layer  salinity  plots,  while  reaching  similar  values  for  the  time  length  of  the 
simulations,  change  at  different  rates.  The  rate  for  simulation  8  is  nearly  linear  for  all 
heat  flux  values  and  is  equal  to  .0013ppt/day  for  equal  to  0.0025  cal/cm^s.  The  rate 
for  simulation  7  is  not  linear  since  the  entrainment  for  this  simulation  is  strong  enough  to 
cause  a  transition  between  freezing  and  melting  conditions.  This  rate  is  equal  to 
O.OOlppt/day  for  equal  to  0.0025  cal/cm  s.  While  this  difference  in  rates  appears 
small,  examination  of  the  plots  shows  that  for  simulation  7  the  rate  appears  to  tend 
toward  equilibrium  and  may,  after  a  long  enough  period,  become  negative  or  decrease  the 
salinity  of  the  mixed  layer.  Meanwhile,  the  rates  for  the  various  heat  fluxes  in  simulation 
8  constantly  increase  the  salinity  of  the  mixed  layer.  Ultimately,  comparing  the  results  of 
simulation  8  with  simulation  7  suggest  that  the  low  ratio  of  OC^gAT  to  P^gAS  from 
simulation  8  cannot  support  as  strong  an  entrainment  rate  as  the  higher  ratio  of  CC^gAT  to 
Pag  AS  from  simulation  7. 
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V.  CONCLUSIONS 


A.  SUMMARY 

The  coupled  ice-mixed  layer  system  is  a  delicate  balance  of  physical  processes 
that  are  dependent  on  both  initial  conditions  and  forcing  conditions.  Changes  in  the  heat 
and  or  buoyancy  fluxes,  wind  speed,  initial  ice  thickness,  and  initial  mixed  layer  depth 
can  all  initiate  feedback  processes  that  can  lead  to  either  ice  growth  or  melting.  This 
balance  of  physical  processes  is  so  tenuous  that  a  moderate  change  in  any  quantity  or 
combination  of  quantities  rarely  produces  an  intuitive  result.  Based  on  this  balance  it  is 
easy  to  understand  why  different  basins  in  the  Arctic  at  constant  latitude  can  have 
different  amounts  of  ice  thickness  at  the  surface.  Similarly,  polynya  formation  and 
maintenance  is  better  understood  in  light  of  these  intricate  relationships.  Ultimately,  ice 
formation  or  melting  cannot  be  forecast  from  surface  forcing  conditions  alone,  as  has 
been  attempted  in  the  past.  Thorough  knowledge  of  the  site-specific  temperature  and 
salinity  profiles  and  fluxes  are  necessary  in  conjunction  with  the  surface  forcing 
conditions  to  predict  entrainment  and  ice  thickness. 

Thermobaricity  is  a  real  and  significant  factor  for  quantifying  the  strength  of 
entrainment  and  the  turbulent  kinetic  energy  of  the  coupled  system.  In  comparing 
simulations  between  different  initial  mixed  layer  depths  the  thermobaricity  seems 
particularly  important  or  sensitive  to  wind  forcing  as  compared  with  changing  heat 
fluxes.  In  other  words,  the  greatest  changes  in  mixed  layer  depths  were  found  for  the 
varying  wind  cases  with  an  initial  mixed  layer  depth  of  1000  m.  For  these  cases  the 
change  in  mixed  layer  depth  could  be  as  large  as  60  m.  The  comparable  case  for  initial 
mixed  layer  depth  of  100  m  could  only  achieve  a  maximum  of  15  m.  (To  see  these 
maximum  possible  depths  compare  Figures  12  and  16.)  This  scenario  is  especially 
important  for  sustaining  polynya  events  and  indicates  that  perhaps  surface  wind  forcing 
does  more  to  maintain  polynya  conditions  than  net  surface  heat  flux.  Furthermore,  this 
mechanism  for  rapid  and  strong  deepening  of  the  mixed  layer  is  important  because  it 
likely  helps  explain  deep-water  regeneration  on  a  global  scale. 


55 


Another  significant  result  of  the  simulations  was  an  observed  relation  between 
greater  ice  thickness  and  higher  heat  storage  in  the  mixed  layer.  In  particular  as 
described  above  every  simulation  shows  that  for  the  greatest  freezing  rate  it  takes  the 
mixed  layer  temperature  longer  to  come  into  equilibrium  and  thus  stores  and  holds  on  to 
more  heat  as  it  comes  into  equilibrium.  Simulation  2  and  simulation  6  (in  Figures  1 1  and 
15,  respectively)  appear  to  be  exceptions  to  the  rule,  however.  For  their  given  conditions 
it  appears  that  the  simulation  did  not  run  long  enough  to  allow  the  mixed  layer 
temperature  to  come  into  equilibrium.  This  can  be  seen  where  the  plots  are  still  rising  up 
and  starting  to  bend  toward  equilibrium.  This  result  is  significant  because  it  indicates 
that  thick  ice  at  the  surface  acts  as  a  cap  on  the  coupled  system.  This  cap  does  not  stop 
entrainment  from  melting  the  ice  at  the  surface;  instead,  melting  begins  before  the  system 
comes  to  thermal  equilibrium  indicating  that  when  there  is  ice  at  the  surface  the  entrained 
heat  works  to  melt  the  ice  before  it  works  to  warm  the  mixed  layer. 

Perhaps  the  most  significant  result  of  the  simulations  was  to  verify  the  complexity 
of  the  interactions  that  occur  in  the  coupled  system  and  which  are  represented  in  the 
equations  developed  in  the  theory  section.  For  example,  as  predicted  by  the  equations,  no 
one  factor  was  determined  to  be  the  most  important  to  quantify  for  determining 
entrainment  strength.  Instead,  because  of  the  feedbacks  inherent  in  the  system,  a 
condition  that  led  to  strong  entrainment  in  one  case  may  not  lead  to  strong  entrainment  in 
another  case,  depending  on  the  other  initial  and  forcing  conditions.  Because  of  these 
feedbacks,  it  is  critical  that  modeling  of  coupled  ice-mixed  layer  system  include  the 
realistic  thermodynamics  described  above.  Similarly,  this  work  has  demonstrated  that  the 
PIPS  model  would  become  significantly  more  useful  for  ice  prediction  if  the  necessary 
physics  were  included. 


B,  RECOMMENDATIONS  FOR  FURTHER  STUDY 

There  are  many  possible  topics  of  further  study  for  high  latitude  mixed  layer 
dynamics  and  thermodynamics.  Preconditioning  for  polynyas  by  advective  processes  at 
depth  like  warm  subsurface  eddies,  as  mentioned  by  both  Hakkinen  (1987)  and  Lemke  et 
al,  (1990),  should  be  explored  as  a  possible  mechanism  to  aid  entrainment  in  polynya 
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initiation  and  maintenance.  From  an  atmospheric  perspective,  further  research  of  the 
effects  of  precipitation  in  ice-free  and  ice-covered  regions  is  important  for  greater 
resolution  of  heat  flux  data  and  the  total  energy  balance.  Net  surface  heat  flux  into  the 
system  needs  to  be  thoroughly  explored  just  as  net  heat  flux  out  of  the  coupled  system 
was  examined  in  this  work.  Similar  sensitivity  studies  of  net  heat  flux  both  into  and  out 
of  the  coupled  system  need  to  be  conducted  for  an  initial  ice  thickness  value  equal  to 
zero.  For  this  study  it  is  critical  that  the  model  be  able  to  transition  freely  and  correctly 
between  the  ice-free  system  and  the  ice-covered  system. 
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